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Abstract 

We study a new quintic discrete nonlinear Schrodinger (QDNLS) equation which 
reduces naturally to an interesting symmetric difference equation of the form n +i + 
4> n -i = F((j) n ). Integrability of the symmetric mapping is checked by singularity 
confinement criteria and growth properties. Some new exact localized solutions for 
integrable cases are presented for certain sets of parameters. Although these exact 
localized solutions represent only a small subset of the large variety of possible 
solutions admitted by the QDNLS equation, those solutions presented here are the 
first example of exact localized solutions of the QDNLS equation. We also find 
chaotic behavior for certain parameters of nonintegrable case. 
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1 Introduction 



Discrete solitons in nonlinear lattices have been the focus of considerable at- 
tention in diverse branches of science [1,2,3]. Discrete solitons are possible in 
several physical settings, such as biological systems [4], atomic chains [5,6], 
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solid state physics [7], electrical lattices [8] and Bose-Einstein condensates [9]. 
Recently, the existence of discrete solitons in photonic structures (in arrays of 
coupled nonlinear optical waveguides [10,11,12,13,14,15,16] and in a nonlinear 
photonic crystal structure [17]) was announced and has attracted considerable 
attention in the scientific community. Photonic crystals, which are artificial 
microstructures having photonic bandgaps, can be used to precisely control 
propagation of optical pulses and beams. They are very useful for optical 
components such as waveguides, couplers, cavities and optical computers. It 
is possible to make discrete waveguides using photonic crystals. In this sit- 
uation, "discrete solitons" appear naturally and have interesting properties. 
Many scientists believe that the discrete solitons can have an important role 
in this technology. 

Here, we consider the physical system described by the discrete nonlinear 
Schrodinger (DNLS) equation. Most relevant for realistic applications (partic- 
ularly to nonlinear optics) is the DNLS equation with a cubic (Kerr) nonlin- 
earity: 

.(h^ + a (^ n+l _ 2^ n + ^, n _ 1 ) + f3\i) n \ 2 i) n = , (1) 



where ip n are complex variables defined for all integer values of the site index 
n. The DNLS equation was used by Christodoulides and Joseph [11] to model 
the propagation of discrete self-trapped beams in an array of weakly-coupled 
nonlinear optical waveguides. In such an array, when low intensity light is 
injected into one, it will couple to more and more waveguides as it propagates, 
thereby broadening its spatial distribution (diffraction). High intensity light 
changes the refractive index of the input waveguides through the Kerr effect 
and decouples them from the rest of the array. Certain light distributions 
propagate while retaining a fixed spatial profile among a limited number of 
waveguides. These are discrete spatial solitons. Experimental results for optical 
waveguide arrays, confirming the validity of the model, have been reported by 
Eisenberg et al. [12]. 

There are many works about the DNLS equation and some its modifications. 
It is well-known that the standard DNLS equation is not completely integrable 
[18]. The integrable discrete nonlinear Schrodinger equation (Ablowitz-Ladik 
(AL) system), 



i ~d~t~ + a ^ n+1 ~ 2 ^ n + ^ n ~ 1 ^ + 7l^n| 2 (^n+l + tpn-l) = , (2) 

was proposed by using inverse scattering method [19]. The AL system has N- 
soliton solutions and a rich mathematical structure. However, it is not phys- 
ically realistic, because it does not contain the Kerr-nonlinearity term. Thus 
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the deformed DNLS equation 



i ~dt L + a ^ n+1 ~ 2 ^ n + ^ n -^ + ^l^™! 2 ^ + F (^n+l, Ipn, 1pn-l) = 0, (3) 

where F(x, y, z) is a polynomial function of x, y, z, was proposed in several 
works. For example, 



F(ip n+1 ,ip n ,tp n ^ 1 ) = 7|^ n | 2 (V>„+i + ^n-l) (4) 

was considered by several authors [20,21,22,23,?]. Equation (3) with Eq.(4) 
can be viewed as a deformation of the DNLS equation (1). Indeed, by using 
the parameterization 7 = (f3' — (3)/2, one can see that, for f3 = (3', Eq.(3) 
with Eq.(4) reduces to the standard DNLS equation (1), while for (3 — > it 
gives the AL system. Kivshar and Salerno explained the physical meaning of 
Eq.(3) with Eq.(4) [21]. They proposed that an application of Eq.(3) with 
Eq. (4) may be found in nonlinear optics in, that the DNLS equation describes 
interactions of partial TE (Transverse electric) modes in array of (focusing or 
defocusing) waveguides. Thus, in spite of the fact that the AL system itself is 
not physical, it can appear as a particular case of a more general and physically 
better justified discrete model Eq.(3) with (4). 

The existence of localized solutions is well known for spatially discrete systems 
(including quintic and higher order nonlinearities) [25,26,27,28]. Moreover the 
stability of these solutions at weak coupling is established as well, independent 
of the integrability of the underlying equations of motion. 

However, nobody has found exact soliton-like solutions of Eq. (3)-type discrete 
systems. Concrete forms of localized solutions are, definitely, useful for ana- 
lyzing real physics. Moreover, from physical point of view, investigating what 
kind of Eq.(3) has exact localized solutions and what kind of solution exits are 
important. Thus the problems we want to address here are: Does Eq.(3) with 
Eq.(4) have soliton-like solutions and what kind of deformed DNLS permits 
soliton-like solutions? 

Here, we propose a new model of the propagation of discrete self-trapped 
beams in an array of weakly-coupled nonlinear optical waveguides, the quintic 
discrete nonlinear Schrodinger (QDNLS) equation 



+7l^n| 2 (VWl + ^n-l) + ^nl^VWl + ^n-l) = 0. (5) 

In physical problems, the quintic nonlinearity can be of equal or even higher 
importance to the cubic one [29] as it is responsible for stability of localized 
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solutions. This equation (in normalized coefficients a — 1,7 = —AS//3) can be 
derived from the Hamiltonian 

H = - J2M* n+1 + «n + i) - y 5 E ln ( X + ^n| 2 + ^l^l 4 )' ( 6 ) 



n 



with the deformed Poisson brackets 

{VV^CJ = *(1 +7l^n| 2 + S\*Pn\ 4 )$nm , {^n, VVn} = {C> = ■ ( 7 ) 

In general, 

rnnl ( OB 8C dCdB\„ . , l2 , (4 , 

{B - C} ^?(«-«)( 1+ *'" |2+ *" 14 '- < 8 > 

The equation of motion is 

j, n = {H,il> n }. (9) 



Surprisingly, the QDNLS equation has exact soliton solutions (bright and 
dark) and Jacobi elliptic function solutions, as we shall show here. 

In this letter, we discuss the integrability of mappings which are reduced from 
the QDNLS equation and exact localized solutions of the QDNLS equation. 



2 Integrability of mappings and exact localized solutions of the 
QDNLS equation 



We begin by looking for solutions of the form ip n (t) = 4> n e~ lu}t where n is 
real. Substitution into Eq.(l), (4) and (5) yields the corresponding symmetric 
difference equations, respectively 



0„+l + n _l 
0„+l + n _l 



-(u/a) + 2)0 n - (P/a)<j> 3 n , a, f3 ^ 0, 
[-uj + 2a)0 n - (3<Pl 



a,/3,7^ 0, 



(10) 



and 



<Pn+l +<Pn-l = 



'-uj + 2a)(j) n - (3<Pl 
a + 702 +5<Pi 



a, 0,6^ 0. 



;i2) 
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Equation (10) was studied by many authors [18] and has chaotic property. We 
apply a singularity confinement (SC) test to these difference equations [30]. 
Equation (11) does not pass the SC test except the case in which Eq.(ll) is 
reduced to a linear mapping n+1 + = {—P/ r i)<j> n ( m the case of (u — 
2a) 7 = a (3). Thus we may conclude that Eq.(ll) is non-integrable except the 
case in which leads to a linear mapping. In the case in which Eq.(ll) leads to a 

linear mapping, it has solutions <p n = v / c -!_ 4 {( C ' + ^ 2 " 4 ) (4>i — c '~ v y 2 ~ 4 0o) _ 
( c ~^ gZi ) W (0i - c+v f^ i o ) } where C = -/3/ 7 and O , <Pi are initial values. 

In Eq.(12), divergences caused by initial values are confined to only one-site 
without losing the initial information for certain sets of coefficients. That 
is, under that constraint, Eq.(12) has the SC property. So Eq.(12) may be 
integrable for certain sets of coefficients. Indeed, this equation has soliton 
type solutions for certain sets of parameters. 

Substituting ip n = ^e~ lujt to Eq.(5), we obtain the multi- linear form 



(u - 2a)/ n+1 /^/ n _ 1 # n + afnf n -ig n +l + a/n+l/nffn-l + Pfn+lfnfn-l9n 
+lfnfn-l9n+l9n + 7/n+l/n#n#n-l + 5fn-l9n+l9n + ^fn+l9n9n-l = 0. 

Using the standard procedure of the Hirota method, we obtain the following 
exact solutions. 

Bright type soliton 

rj) n (t) = Asech(nlog(p) + log(n )) e -^, (13) 



where u Y & ' P ~ 2^ ' A ~ ± 2y aspa-w) ancl n ° 

is arbitrary constant. 

Dark type soliton 

if> n (t) = Atanh(nlog(p) + log(n ))e" ia; ', (14) 

where u = ahAEgli^ , p = y iS5g5ir , A = ±V=I5 and „„ is 
arbitrary constant. 

The QDNLS equation also has Jacobi elliptic function solutions for certain 
sets of parameters. 
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Fig. 1. Bright soliton of Eq.(5). a = 1,(3 = 0.2, 7 = 4.78701,(5 = -0.2. Time in this 
snapshot is 100. 

Here, we show Fig.l, which is the result of numerical evolution of an initial 
exact bright soliton. We should note that the above localized solution is stable 
for certain parameters. 

What happens for the sets of parameters for which Eq.(12) does not pass the 
SC test? We show two figures ( (<f> n+ i,<f> n ) phase-plane plot) here. As can be 
seen, Fig. 3 uses a set of non-integrable parameters and chaotic behavior is 
obtained (compare to Fig. 2). 

Let us explain why the mapping (12) is integrable for certain sets of pa- 
rameters. It is well-known that the McMillan mapping [31] z n+ \ + z n -\ = 
(A + Bz n )/(C + Dz%) has a conserved quantity and could be solved by Ja- 
cobi elliptic functions [32]. Equation (12) for certain sets of parameters re- 
duces to the McMillan mapping by cancellation of common factors in de- 
nominator and numerator of Eq.(12). These cases of certain sets of param- 
eters are equivalent to the cases in which Eq.(12) passes the SC test. For 
example, in the set of parameters of Fig.2, Eq.(12) reduces to the McMil- 
lan mapping n+1 + <p n -i = — 70 n /(— 22 + 110^), which has an invariant 
—22(0^ + 0^_ x ) + 70 n 0„_i + 110^0^_ x and Jacobi elliptic function solutions 
[32]. 

We note that the quintic NLS equation with fourth order dispersion iu t +au xx + 
tbuxxxx + c\u\ 2 u + d\u\ 4 u = where e is a small parameter, corresponds to a 
continuous analogue of Eq.(5) [33,34,35]. This equation has exact bright (i.e., 
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2. A collection of orbits of the map generated by Eq. (12).(Integrable case: 
37, a = 22,(5 = 7, 7 = —33,5 = 11. These parameters were determined by the 
test. Each orbits are described by Jacobi elliptic functions.) 




3. A collection of orbits of the map generated by Eq.(12). (Non-integrable case: 
37, a = 22, (3 = 6, 7 = -33, 8 = 11.) 
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A sech(kx)e lt ) and dark (i.e., Atanh(/cx)e rf ) soliton solutions for certain sets 
of parameters [34,35,36]. The existence of exact localized solutions originates 
from the balance between quintic nonlinear term and 4th-order derivative 
term. This equation appears in the optical fiber transmission studies [34,35,37]. 

We can easily generalize to the following DNLS equation: 



+72|^n| 2 0n+l + V'n-l) + Aj^nlVn 
+74|^n| 4 0n+l + 1pn-l) + At^nlVn 
+ ••• 

+l2N-2\^n\ 2N ~ 2 {^n+l + V'n-l) + l^nl^" Vn 

+72^v|^n| 27V (^n + l+^-l) =0. (15) 

Using a solution ansatz of the form ip n (t) = (p n e~ luJt , where <p n is real, sub- 
stitution into equations (15) yields the corresponding symmetric difference 
equation 



, , . _ y-uj + 2a)<p n - fofi - P2N-2^ n N - 1 nfi , 

n,,+s n "- 1 ' «+ 72 ^+ 7 ^+.--+ 7 ^_ 2 0r- 2 +72^r • ( ' 



This mapping has the SC property and low-order degree growth property 
when the coefficients satisfy some constraints. In fact, it is also possible to 
construct some exact solutions by using standard Hirota method. We may 
consider the existence of these exact localized solutions also originates from 
the balance between 2N+1 order nonlinear term and 2iV-order derivative term 
which comes from the higher order term of Taylor expansion of second-order 
difference term. 



3 Non-integrable mappings and their growth properties 



Next, we consider the forced quintic discrete Nonlinear Schrodinger (fQDNLS) 
equation 



+5\^ n \\^j n+1 + V>„-i) + / exp(-iwf) = 0, (17) 
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Fig. 4. A collection of orbits of the map generated by Eq.(19). 

where / is a constant parameter. This equation can be reduced by ifi n = 
(p n exp(—iujt) into 

i, + <t, _ / + (-■" + 2 «)^- - Wji no, 



Our problem is whether this is integrable or not. After we apply the SC test, 
we note Eq.(18), for certain sets of parameters passes the SC test without 
losing initial information. For example, we can consider the simple mapping 

n+1 + n _ 1 = 1/(1-0*). (19) 



This mapping has the SC property. But this result is a counterexample to 
the physical intuition that nonlinear equations containing forcing terms are 
chaotic. To verify this, we perform numerical calculation on Eq.(19). The result 
obtained from Eq.(19) is shown respectively as ((fi n +i, (fin) phase-plane plots in 
Fig. 4 for varying choices of 0o,0i- (In contrast, (0„ + i,0 n ) phase-plane plots 
(Fig. 2) of Eq.(12) show periodic orbits.) 

This situation is similar to the case of Hietarinta and Viallet [38]. Algebraic 
entropy is useful to check integrability in this case [39,40,41]. The main ar- 
gument is that a generic, non-integrable mapping has an exponential degree 
growth, while integrability is associated with low growth and is typically poly- 
nomial. The degrees of the iterates of Eq.(19) are 0,1,4,17,64,241,900,..., i.e. a 
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exponential degree growth. Thus we can assert that Eq.(19) is non-integrable. 
This fact can be confirmed by considering cancellation mechanism of common 
factors in denominator and numerator of Eq.(19) (or Eq.(18)). The denomina- 
tor and numerator of Eq.(19) (or Eq.(18)) do not have common factors. Thus 
Eq.(19) (or Eq.(18)) cannot be reduced to the McMillan mapping. 

From above analysis, we can easily produce the following generalized mapping: 



a + ai(f) n H h a N ^i<j)^ 1 

b + h<j) n + 6 2 02 + ■ ■ ■ + b N <f$ 



0n+i + (f>n-i = i I —2 I n— In- ( 2 °) 



This mapping also passes the SC test for certain sets of parameters, but has 
exponential degree growth. Thus (0 n+ i,0 n ) phase-plane plots of Eq.(20) also 
trace chaotic orbits. 



4 Conclusions 



We have analyzed the QDNLS equation. The QDNLS equation is not com- 
pletely integrable, however the QDNLS equation has a number of exact lo- 
calized solutions exists to the QDNLS equation provided that coefficients are 
bound by special relations (i.e. partially integrable). The set of localized solu- 
tions include particular types of solitary wave solutions, dark soliton solutions 
and periodic solutions in terms of Jacobi elliptic functions. Clearly these local- 
ized solutions represent only a small subset of large variety of possible solutions 
admitted by the QDNLS equation. Nevertheless, the solutions presented here 
are found for the first time and they might serve as seeding solutions for a 
wider class of localised structures which exist in this system. We hope that 
they will be useful in further perturbative and numerical analysis of various 
solutions to the QDNLS equation. 

We mention traveling wave solutions of the QDNLS equation. Flach et al. stud- 
ied moving pulses of the generalized discrete nonlinear Schrodinger equation [42] 



+G(|^„|> n = 0. (21) 

Their result shows that the AL system is only generalized discrete nonlinear 
Schrodinger equation with sech-type moving pulses, i.e. the QDNLS equation 
does not have exact sech-type exact traveling wave solutions. Although the 
QDNLS equation does not have exact traveling wave solutions, there is still the 
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possibility that non-exact traveling wave solutions exist. Such solutions must 
conserve the integrals of motion (number of soliton, energy and momentum). 

We have also found interesting properties of mappings derived from the de- 
formed DNLS equation. 

Further studies (e.g. stability of exact localized solutions and bound states of 
multiple pulses, and the existence and properties of travelling wave solutions) 
will be addressed elsewhere. We believe our result can be applied to the study 
of discrete optical (and other) solitons. 

KM thanks Y. Kivshar, N. Akhmediev, A. Ankiewicz, K. Kajiwara, S-I. Gotou, 
M. Oikawa and T. Tsuchida for useful suggestions. KM was supported by a 
JSPS Fellowship for Young Scientists. KM and YO thank University of Sydney 
for warm hospitality. NJ's research was supported by the Australian Research 
Council. 
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